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ABSTRACT 


The paper is concered with a theoretical study of shearing flow 
bounded by a wavy surface with consideration of turbulent flow above 
the air-sea interface. Account is to be taken of turbulence in the flow 
through the use of Reynolds stresses associated with turbulent flow. 

An adaptation of the mixing length theory as applied to pipe flow is made 
for channel flow and the resulting mixing length versus energy relation- 
Ship is incorporated in the Reynolds stress term. Finally, the rate of 
wave growth is calculated from the normal surface pressure in phase 
with the wave slope. 

A curvilinear coordinate system which follows with the wave train 
is used in order to simplify the formulation of the problem. All para- 
meters are non-dimensionalized and the analysis is made considering 


a velocity profile adapted from pipe to channel flow. 
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SYMBOLS AND ABBREVIATIONS 


a non-dimensional wave amplitude 

A constant parameter 

e non-dimensional wave celerity 

C constant 

d position at which upper boundary conditions are applied 
I turbulent energy 

E wave energy 

E rate of energy input to wave 

F function determined by solution 

i Froude number 

g acceleration of gravity (980 cm sec) 

Hi channel height 

J Jacobian 

K kinematic eddy viscosity 

p mixing length 

oY mixing length (perturbation term) 

L mixing length (zero-order term) 

| pressure 

q velocity of main flow above boundary layer 
U velocity of main flow in the channel 

u non-dimensional velocity parallel to 5 axis 
u,. friction velocity 


< 


Se ee ee Im ER 


non-dimensional velocity parallel to ‘| axis 


COnStane 

real component of pressure 
imaginary Component of pressure 
real component of shear stress 
imaginary Component of shear stress 
constant parameter 

mixing length constant (0.4) 
non-dimensional wave length 
stream function 

molecular viscosity 

air density (0.0013 g cm73) 
water density (1.027 g cm73) 
Shear stress 


curvilinear coordinates 
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1. INTRODUCTION 


A semi-empirical study of turbulent air flow over a wavy surface 
is made in which normal and tangential stresses on the surface are 
examined with the object of estimating the growth of water waves due 
to a turbulent wind. 

Ursell (1956) gives a summary and a critical review of the theories 
of wind-wave formation extant in 1956. The review gave the necessary 
impetus for further investigations, both theoretical and experimental, 
in this area. Since 1956 there has been much work done on the generation 
of water waves by surface stresses induced by the wind. Miles (1957) 
and Benjamin (1959) considered laminar air flow over a wavy water 
surface and deduced formulae for exponential wave growth due to an 
instability mechanism. Later Phillips (1966) included air-stream 
turbulence in the instability model from a semi-empirical point of view. 

In his initial work, Miles (1957) assumes an inviscid, incompress- 
ible, two-dimensional, parallel shear flow coupled with a prescribed 
two-dimensional deep-water gravity wave. Turbulent fluctuations in 
the perturbation equations of motion are neglected, even though they 
are decisive in maintaining the mean shear flow. Furthermore, the wave- 
induced velocities and pressure are small, which justifies linearization 
of the equations. The water is also assumed to be inviscid, incompress- 
ible, and irrotational, and the slope of the displaced surface is assumed 


to be small. Only that component of the aerodynamic force in phase 
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with the wave slope is important in the wave-generation mechanism. 
The viscous drag forces of the air are assumed to be negligible compared 
with the normal pressure in their effect on the surface wave. 

This model is used to determine a first-order approximation to the 
disturbed motion of the air and the consequent energy transfer to the 
wave. The velocity distribution is taken to be that for turbulent flow 
Over an aerodynamically smooth water surface. The result rests on the 
solution of the inviscid Orr-Sommerfeld equation which has a singularity 
at the height where As = C (the critical level). An exponential 
increase in amplitude of water waves with phase velocity c occurs if 
the mean air flow velocity Uw in the direction of wave travel varies 

1 
with height above the surface such that Ws at the critical level is 
negative. This is almost always the case. 

Phillips (1966) attempts to include wave-induced turbulence in the 
instability model by including in the momentum flux equation an 
appropriate integral over the entire layer of air of the contribution due 
to perturbation eddy viscosities. This integral is multiplied by an 
unknown constant A which is estimated from experimental results from 
flow over fixed wavy surfaces. The uncertainty in A is estimated to 
be + 50 percent. Phillips concludes that even though there exists 
a high degree of uncertainty in A the results do provide a framework 
in which experimental measurements may be interpreted. 

According to Miles (1967), experimental results indicate that a 


theoretical model based on laminar flow may be adequate in the 
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laboratory but not on an oceanographic scale. The inviscid, laminar model 
underestimates the energy transfer from wind to waves over at least a 
significant portion of the spectrum for an open sea. This suggests that 
the significance of wave~-induced perturbations in the turbulent Reynolds 
stresses for momentum and energy transfer from wind to wave must 
increase with an appropriate scale factor. Miles concludes that further 
theoretical work in this area appears to demand some ‘ad hoc’ hypothesis 
to describe these perturbation Reynolds stresses. 

This can be done by considering a slightly generalized form of the 
classical mixing-length hypothesis. A means must be found to predict 
the value for the change in mixing length upon which the perturbation 
Reynolds stresses are based. Zagustin (1963) has given a semi-empirical 
solution for turbulent flow in a pipe which gives good agreement with 
experimental data and which predicts the value of the mixing length. 

This model is generalized below to include flow over a wavy surface in 

a manner analogous to that of Miles and Benjamin in that it assumes 

the flow to follow the wavy surface. Viscous stress terms are replaced 

by the Reynolds stresses associated with turbulent flow. Normal Reynolds 
stresses are neglected. Zagustin's method of closing the system of mean 
turbulent equations of motion is not unique (see, for example, Squire, 
1959), but there seems to be no good reason to use any other model in 

its place. The investigation below details the application of Zagustin's 


work to this model and outlines the solution of the differential equations 


FS 


using numerical methods and appropriate boundary conditions. The 
unknown parameters in the model are fitted to experimental data describing 


the normal pressure in phase with the wave slope. 
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2. PORMULATION OF THE MODEL 


The zero-order flow is taken as flow in a two-dimensional channel. 
All lengths are made non-dimensional using the channel halfwidth, H. 


Velocities are made non-dimensional using the midstream velocity, U, 


is 
All stresses are made non-dimensional by dividing by Ui , where 
Io. So. 


is the density of air, 
2.1 THEZERC-ORBER FLOW 


The velocity is expressed in terms of a stationary and a fluctuating 


(turbulent) part: 

(2.1) — Sa. 

Following Zagustin, the Reynolds shear stresses T are written as 
(2.2) ct = CS qradq) 


where q is the mean speed. This can be considered a defining 
equation for L (although here it appears as the usual mixing length) 
in that I will not be prescribed, but will be part of the solution. Then 


the two-dimensional mean equation of motion is 
N 
(2. 3) (U e grad) U = -gradP + Ik x grad iq 


A 
where IK is a unit vector normal to the plane of the flow, P is the 


2. 
pressure, and K = , grad q is the eddy viscosity. The turbulent 


ee 


kinetic energy, E' = Wee wu , 1S assumed to have the form 
(2.4) ial | 2 (grad q)2 


where © is an unknown constant. It is assumed that the rate of 
dissipation of E' is proportional to E'. 

A flux of turbulent kinetic energy into any volume element is assumed 
to occur when the quantity L varies over the element. Thus; tite 


rate of inflow of E' is 
(25) x| | (ora |" / n) ds 
S 


N 
where n is the normal to any surface and © is an unknown constant 
depending on the turbulent intensity. The turbulent kinetic energy is 


constant at any point, and the above assumptions then imply that 
(2 . 6) \ { (rad i> acc °})j E' dv = 0, 
S 


where B is an unknown constant. Applying the divergence theorem, 
considering V to be an arbitrary volume, and using equation (2.4), 


equation (2.6) becomes 


(2.7) w2|. +8 L2radq)? = 0. 


Equations (2.3) and (2.7), together with the continuity equation, 
form a closed set for the mean velocities and the quantity E . The 


constant B must be determined experimentally. 
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For flow between two parallel walls as shown in Figure 1, the 


continuity equation is trivial, and equations (2.3) and (2.7) become 


(2. 8) -grad P + d_ oe id U | E10 
dy dy 
(2.9) ais 8 | “a yh ye = 0. 
dy2 dy 





Figure l. 


Defining B =A for the zero-order flow case and setting 
Z 


equation (2.9) becomes 


(2.10) a2 | n, a eg 
and equation (2.9), after integrating, becomes 
(2.11) hy = | “@ Uy 
dy 
Equations (2.10) and (2.11) give 
(Ga) a2 L 


A 
= -Ahy, 
dy 2 2 = 


lbs, 


and 


, =| grad P|, 


which has the boundary conditions: 


al ns 


y= 0 


- 


where | o is a quantity analogous to the roughness length in mixing- 





y—-4 


length theory. Integrating equation (2.12) twice and applying these 


boundary conditions results in 


(2.13) | = 2D (1 - y3) 4 = 


Next, selving equation (2.4) 0 ferd U and substituting for iz 


dy 
gives 
(2.14) d U ee Jhy ae 
Ms Ah 3 
ei (] - 
ony 2 | v 


By defining 


pe Pe 2 and 


equation (2.14) simplifies to 





=| 
(2.15) 12 alll a va , dz. 
A Jhb a2 
ae S 
Integrating equation (2.15), applying the no-slip condition atz= - Vb; 


and returning to the y coordinate gives 
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8 | all =) 3/2 


The mean velocity at mid-channel is 


WU (y) 





ye= 0 
which, when applied to equation (2.16), results in the expression 


(2.17) —— ie 


(2.18) Maga = 1 - ae 


The expression (2.13) for a can be simplified by recalling that 
von Karman's constant, 3 , gives the rate of growth of | at the 


wall. Therefore, 


(2.19) UL = al 








where K ~~ 0.4 is the accepted value. Then, from equations (2.13) 


it is found that 


(2.20) Ah = 1.6 

(2.21) Bb = Slee 735 = 
S 

(272) & Se a . 


Following Benjamin (1959), the equations of motion in terms of E 


and ng are 
(2. 30a) (bs, ag = L, ee oe 1/2, 4s + G42) = - 


‘ 
b) 1b eg . 4 bh ve 1/24 (4? 42 ye - il +b, 


where the subscripts denote partial differentiation. 
The stream function is composed of a zero-order term plus a 


periodic perturbation. When there is no wave, the stream function is 


y 
(Zech) V4 = | [Uoo < | a 


where ==" y —ole 
iv (y) is the velocity profile from equation (2.18) expressed in 


curvilinear coordinates 


tanh! | u an, hdl 
(2. 32) T= ee el 


When the water surface is disturbed by a stationary wave, the stream 


function becomes 


(2.33) Ww (2) = We F(W ) + LU ee jam oikg 


where F(Y) is to be determined. 


The mixing length is expressed in a form similar to that for 7 ; 


(2.34) dX = Z (y) + a Yin elt 
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where | (1K) is given above and YX (¥ ) is to be determined. 


Similarly, 


(2.39) B = o + ad alks 


where A is given above and oo is an unknown parameter to be 
determined empirically. The velocity components parallel to 3 and ” 


are given by 


(2.36a) u = yyy = U -cta Gace + U ire} eik§ 
4 -ika { em + LUM = a || et ave 


To formulate the differential equation in Y. to be solved, 


(oy 
— 
i 

3 
iI 


substitution is made for YL , B and qin ~quation (2.7) expressed in 


curvilinear coordinates: 
(267) 7 WV ee =o 
yi 
Simplification by elimination of all terms of O(a) and above yields 


(2738) P f JX @e 4) + AHF" + &D + 


Neu (Eo koh = 6 


where 


Zo 


Similarly, when equations (2.30) are expanded, cross differentiated 


to eliminate pressure terms, and linearized in a, one obtains 
(2.39) ie _ i ee ellie a 
5 ( U ey Uae le) 0) F sgn( [J ) = 
= k*) (HE" + p Y ) + 
ui 
e “Yl ig" - 2kG' + 2k%G) 


where 


?) 

q 
4 
TS 


a= UL 

and where Lule is positive for the lower half of the channel. Note that 
a singularity occurs at mid-channel ( UD = 0). Equations (2.38) and 
(2.39) are to be solved for F and We ‘ 

The boundary conditions are as follows. At the interface where 
vy = Q, U (0) = 0 to satisfy the no-slip condition, and v = 0 because 
the wave is stationary. Therefore, from equation (2.36b): 
(2.40) F (0) meSeas 
The velocities in the water must take on negative values in order that 
the shearing stress across the interface be continuous. The variable 
component of the tangential water velocity associated with an irrotational 


wave is of the form 


(2.41) wease (SLs 
ein 
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The total tangential water velocity is 


(2°42) u= -c - a Jk’ ihe elk § ; 


Defining the Froude number: 





2 
(2.43) a = Us , and invoking the no-slip condition, the 


tangential air velocity is 





(2.44) Us “eer = ~.2 r oik§ at " = 0, or 
(2.45) ie. -C 2B alae at | = 0 
where B = — . From equation (2.36a) the second lower boundary 


condition then becomes 


(2.46) F'(0.) = - UF (0) + p 


The perturbation in DL at the water surface, L. (0), must also be 
fixed. This must be regarded (along with o. ) as an unknown parameter 
in the problem, to be adjusted empirically. 


The perturbed flow must die away at some large distance (relative 


to k~+) from the wavy surface. Thus, as 4 becomes large, 
(2.47) Fy) ——-— BG 5. *% 

(2.48) F'(y )——> ay) ok 

(2.49) y (y)—>- 0. 
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2.3 NORMAL PRESSURE, SHEAR STRESS AND WAVE GROWTH 


The pressure on the wavy surface is the sum of that due to normal 
atmospheric turbulence and that generated by air moving over the 
undulating surface. The latter pressure predominates in the generation 


mechanism considered here and is expressed as 
(2.50) P(0) = a(%& + idg )elk& 


where oe and & are to be determined. Differentiating equation 
(2.50) with respect to g , substituting the result into equation (2.30a), 


and simplifying yields 
| _ ais z 1 . 2 Ld J 
(2.51) P(0) = Ut CMU eye aa (HF 


+ EY ) + (G' - oe) | 
V¥is=0 


which gives “” and s , 


In the same manner, the shear stress at the wall is expressed as 


(2.52) Tow = a(€ +i € )eik® 


where € and € are to be determined. To evaluate the shear 


stress, equation (2.2) is expanded in terms of Wy 


(2.53) T = 220? Wy a + In YY, Yan) 


SUootituting fOr X ar Ly , or their derivatives where required 


obtains 
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(2.54) Tv @ = | 2 LULU +7 ] 
y = 
which gives € and € . The shear stress seems to be relatively 
unimportant to wave growth, however, and will not be considered herein. 
To compute the wave growth rate the normal pressure in phase with 
the wave slope is required. Expansion of equation (2.50) and comparison 
with equation (2.25) for the wave form show that a positive A results 


in wave growth. The average energy input due to normal pressure is 


» 


(2.55) paces & pd 
mn 


: 


or, in dimensional terms, 


(2.56) y 2 oe U,? S, 


Q 


2 


where W = gk is the wave frequency. 


The surface energy of the wave is 
(2.57) E =1/2 Fur ga? 
where Gus is the water density, so that 
(2.58) a = for gada 


Equating equations (2.56) and (2.58) and solving for oa gives 


2 
(2.59) s3 a "i Ue" $y 


ao 


Integrating produces the familiar exponential wave growth equation: 


(2. 60) a= Saree | = c Ga t 
us 


The rate of energy input due to surface shear-is 


7 » 

(2.61) E=j2 FC. as 
— 2 

(2. 62) EB sae ® eae, 


This energy is expected to go, at least initially, into a surface shear 


flow. 
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3. SOLUTION OF THE FIRST-ORDER PROBLEM 


A solution for the function F and its derivatives must now be 
found so that values for normal pressure, shear stress and wave growth 
rate may be determined. There are several numerical integration methods 
available for this purpose, having the common property of starting at 
one boundary and ending at the other. However, the problem as 
formulated has conditions at both boundaries, and a method of super- 
position of solutions must be used. 

Equations (2.38) and (2.39) can be rewritten as six linear first- 


Order ordinary differential equations. That is letting: 


(3. 1a) F =Y 


} 
b) FE =Y¥, 
Cc) F" = Y., 
d) ee 


gives six ordinary differential equations: 


(3 : 2a) Yy, = Yo 


I 
b a = 
) 5 3 
Cc) Y = 
3 4 


il 


| "N 
d) ee 44 4[ 0 Cl ha = iwi | = 2H'Y, 


v\ 
- (a + kay, - EX, -2Ey, - (e" +kEVy, 


eis WGce 17 on } 


and 


ee ae 


There are now six equations, each requiring an initial input at the lower 
channel boundary. Linearly independent solutions are found by starting 
repeatedly at the lower boundary with different initial values, Y(0). The 
final solution is obtained by superposing the independent solutions to 
fit the required conditions at the upper boundary. 

The independent solutions, denoted by the Superscript j, have 
initial values given by | 


(3. 3) Yp (0) = 4 ! sat sel aen 


1 


where | is the Kronecker delta. The solution to the homogeneous 
form of equations (3.2) can be written as a sum of independent solutions 


to the homogeneous equation: 
4 = > _y 
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where the J's are complex constants. Since 

F(0) = F'(0) = WY, (0) = 0 
for the homcgeneous system, the constants Jy. Jo, and J. in equation 
(3.4) must equal zero. The complete solution is achieved by adding to 
the homogeneous solution a particular solution of equation (3.2). For 


this particular solution, the initial conditions are 


(3.5a) yP = ¢ 


1 
pe E 
oe vs G U 
c) yP = 0 
3 
d) y" = © 
4 
e) - = Xo) 
f) Yo —=e 
6 


= 3 6 
(3.6) -¥(M) = 1g ¥8Y) + VOT) + LQ) + KD), 


where the J's are to be determined. 

The position of the upper boundary, where conditions (2.47), (2.48), 
and (2.49) hold, is not known, and must be determined by a convergence 
criterion. The solution is determined by applying the upper boundary 
condition at various increasing values of 7 = d, solving the resulting 
equations (3.6) evaluated at d forthe J's, and evaluating the wall 


pressure each time using equation (2.51). When the wall pressure 


determined in this manner converges, the procedure is stopped. 
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4. THE,COMPUTE RAS OLUTLON 


The solution for the function F and its derivatives has been 
obtained by solving the problem as outlined in section 3 on an IBM 360- 
67 computer using complex double-precision arithmetic. Major sub- 
routines are incorporated in the program to calculate the velocity profile, 
perform the integration using the Runge-Kutta Gill fourth-order method, 
compute the shear stress, calculate the wave growth rate and evaluate 
the stream function. Appendix I contains the format of the program and 
a set of results obtained. 

The accuracy of the Runge-Kutta Gill integration method depends 
on the size of the step used in the iteration. The value i0°° was 
selected in order to determine whether the surface pressure converged 
on a particular value. Because surface pressure did not converge for 
this step size, it was decreased until convergence was achieved. A 
step size of 10°° was necessary to obtain convergence of the wall 
pressure, and the program running time was about 112 minutes. A final 
optimum step size of ie LOL 4 less than 0.3 and 2.5 X Lome for 4 
greater than 0.3 was then adopted. This step size gives convergence 
on a result essentially the same as that for 107° and a running time 


of 65 minutes, 
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5. RESULTS AND CONCLUSIONS 


The zero-order velocity profile is given by equation (2.18) which 
has ly (0) as a parameter. A plot of U (y) for several values of L (0) 
appears in Figure 3. From the figure it can be seen that as the parameter 
7 (0) decreases, the shape of the profile changes from nearly linear to 
almost a step function. A representative logarithmic profile is also 


shown in the same figure so that the two profile shapes may be easily 








compared.: 
1.0 % 
| Inverse Hyperbolic Tangent | 
and 
Logarithmic Profiles 
0.8 


i (0) = 107 eee 
| () = 107%, 
0.6 aie TURE 


L =o. 


Typical logarithmic 


Joi profile; 





Figure 3, 


0 


The remainder of the analysis is concerned with the determination 
of suitable values for the parameters | (On YL (0), and x . FoR 
the channel under consideration: 

H = 10 meters 

U= 10 meters/second 
which results in a Froude number of 1.018. For the wave number; 

k = 15, 
the instability mechanism considered is important. These values for 
the Froude number and the wave number have been used throughout the 
analysis. 

The following parameter values were first used: 


|.) Dy PDS IM 


Dv C= 
oY = 10. 


Z 


The wall pressure failed to converge for these values. Convergence 
was assumed to have resulted if wall pressure values were nearly 
constant over a significant range of d (say 0.1) or reached a minimum 
value at some point in the channel. (This latter criterion is not 
completely defensible, but the singularity in equation (3.2d) at y = 1, 
itself a peculiarity of the Zagustin assumptions, is probably responsible 
for the erratic behavior of the wall pressure as 4->~1 ~) o~ was then 
varied from -1000 through zero to +1000 and given imaginary values as 


well. The wall pressure changed only in the fourth or fifth decimal 
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which was considered insignificant, although QL changed markedly. 
Thus, the solution seems insensitive to e 


New values for the parameters were then used: 


| (0) = 5¥ 10> 
Lo) = 0 
= pie 


The wall pressure converged for these values, but the value was too 
high in magnitude and of the wrong sign. A value nearer to 5 (Miles, 


1957) was considered appropriate. 


Lo 
Y, (0) 
a Soe 


the resulting wall pressure was positive but still too high in magnitude. 


With 


II 
on 
~ 
— 
= 

t 


| 
= 


Figure 4 summarizes the results of varying la (0) and ma (Oe 


d (position where upper boundary 
1.0 conditions are applied) 






Lo = 2.22X 10° Biot | (0)=5xX10” 


Xo) 





I 


Il 


= 10 


S 101 
Wall Pressure, 
Figure 4, 


of Ge 


It appeared that the wall pressure was more directly dependent on 
y. (0) than on any other parameter. Subsequent investigation showed 
that for any yi the wall pressure was a linear function of the real or 
imaginary part of WL (0) providing that the other component was held 


constant. For 


S 


Lo) 
X 0) 
Ss = 10 


the resulting minimum wall pressure was 3.83 at d=0.48. Figure 5 


5X10 


0.65788 


shows the resulting wall pressure in this case versus d. 





Tae 

a ES oe ree 
0.8 

Vo) =Qves7 eeu: 
ae 

10 

0.4 
0,2 


Wall Pressure, A, 


Figure =o. 


X* was then given a value: 


s* = 1000. 
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No significant change occurred in the wall pressure. Convergence, as 
defined above, always occurred near d= 0.5. 
Table 1 summarizes the values of the wall pressure, shear stress, 


and wave growth rate versus Y. (0) at d= 0.5 for parameter values: 


lo) = 5§x10° 
x = 10. 
Out-of-phase In-phase Wave Growth 
Y. (0) Wall Pressure ,4 Shear Stress,€ Rate, Sa. E S$ 
(aJ 
0 OL -6495.3 =~ 20m. Aas = Eo 
ez a OL 4 oil Ore -14.89 = Mae 
O24 +. 0i -2543.6 Omow —760 
imo + Oi =O ne 14.87 -0.141 
0.8 + Oi 1408.2 PAS| IE 0.349 
iO + 01 3384.0 44.64 0.838 
0.65788+t0i 3.83 ee eal, 0.001 
0.0 + 0.21 -6495.3 = eG = Leu 
O70; + 0.41 =6401-¢ Bee (ae -l 98 
UAO a 0 6l wee aoc Us a3 al She Uae: Lass 
GaGa Oscl @-6214 56 pe /S)5 (he: ee 8/2 
Oe te Wn Ol O27 lees —2 eee — lad 
O265 76610. 11 Oo S —29eee Oe01s 
Table l. Wall pressure, shear stress, and wave growth rate 


as a function of the perturbation surface roughness. 
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It can be seen that imaginary values for YY, (0) have no effect on the 
Shear stress. The required value for LY, (0) is 0.65788 + Oi for a wall 
pressure of 3.83, but the shear stress seems high. (However, the shear 
stress value is a first-order contribution to the total shear stress and is 
reduced by the fractional value of the wave amplitude.) Exponential 
wave growth rates, due to the wall pressure alone, appear reasonable. 
To reduce the shear stress to a more acceptable value in relation to the 


normal pressure, parameter values 


Lo) 
DY. (0) 
y+ = 10 


were used, and give at = ee 


5x 107° 


Ome = 54554 


Normal pressure = 5.0 
Shear stress = 0.001 
Growth rate —=— 0-001. 
However, the wall pressure failed to converge, and had the form shown 
in Figure 6. 
Equation (2.33) for the stream function has been evaluated with 
the values of F calculated from the constants J determined at d = 0.48 
where the wall pressure was found to converge. Parameter values 


were 


5 


Lo 
¥ (0) 
Ss = 10, 


ae 10m 
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Figure 7 shows a plot of the streamlines with the 'cats eye' present 

and centered about 30 degrees ahead of the wave crest. The ‘cats eye’ 
or closed loop is indicative of the wave induced perturbations in the 
flow and characterizes the region of turbulence (Phillips, 1966). Higher 
in the channel the streamlines become less wavy and approach straight 
lines. Table 2 shows values of the stream function near the critical 
layer. 

In view of the foregoing limited results, the following tentative 
conclusions can be made: 

1) The inverse hyperbolic tangent velocity profile, a result of 
the theory of the model, can be fitted to almost any known profile type. 
The profile which has a large value for —U . near the critical level 
and for which the parameter | (0) has the value 5X 107° produces 
surface pressures which can be made to converge on the desired value 
by varying the wall perturbation in Dy . 

2) Wall pressure is extremely sensitive to very small changes in 
2X (0) and insensitive to moderate changes in in ; 

3) Reasonable wall pressure (and wave growth rates) can be obtained 
with an appropriate choice of LV. (0). This indicates that small, unequal 
variations in the surface roughness, such as would be caused by 
parasitic capillaries, are very important in the wave-generation process. 

4) Shear is independent of imaginary values of as (0). Therefore, 


to vary the shear stress, (Ro ( ay (On) oor | (0) must be varied 


individually or in combination as indicated by further analyis of the model. 
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Based on the above results, the following suggestions are made: 
1) | (0) should be increased slightly to 
|) =e 
in order to reduce the shear stress relative to the wall pressure. One 
computer run indicates that the shear stress is halved for this value. 

2) Additional wave numbers should be investigated, say, in the 
range5 < k € 20, and their effect on wall pressure and shear stress 
noted. 

3) The model should be checked using data from Stanton and 
Motzfeld as quoted in Ursel? (1956), and with recently obtained laboratory 
data (see, for example, Zagustin, et al., 1966). 

In summary, the model appears to be fairly well suited toa 
theoretical study of air-stream turbulence and wave generation. 
Unfortunately, the long time required for a computer run prevented the 


calculation of further data using appropriate parametric values. 
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0.025 OQ 312.5 
Om0 Se 0.424 
0.045 0.506 
Q-a0 
5 0.0 
0.002 sys 
0.005 0 W457 
0.008 0.456 
0.10 pas BORE 
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0.020 Omer 
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U035 0.8 1u7 
0.045 Lee 2ele2 
Table 2. 
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